# Forelesning 10, fredag 10.02.2017

#2# Medisineffekt?
x<-matrix(c(23,7,18,13),2,2,byrow=T)
chisq.test(x)$exp
chisq.test(x)

#3# Hardy-Weinberg likevekt?
X <- c(51,42,7)
n <- sum(X)
phat <- (2*X[1]+X[2])/(2*n) #allele frekvens
phat
Phat <- c(phat^2,2*phat*(1-phat),(1-phat)^2) #Genotypefrekvens
Phat
n*Phat   #Forventningsverdier
Q <- sum((X-n*Phat)^2/(n*Phat))
Q        #Ser det ut som forkastning av H0?
pchisq(Q,df=1,lower.tail=F)